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We consider globally regular and black hole solutions in SU(2) Einstein- Yang-Mills-Higgs the- 
ory, coupled to a dilaton field. The basic solutions represent magnetic monopoles, monopole- 
antimonopole systems or black holes with monopole or dipole hair. When the globally regular 
solutions carry additionally electric charge, an angular momentum density results, except in the 
simplest spherically symmetric case. We evaluate the global charges of the solutions and their ef- 
fective action, and analyze their dependence on the gravitational coupling strength. We show, that 
in the presence of a dilaton field, the black hole solutions satisfy a generalized Smarr type mass 
formula. 
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I. INTRODUCTION 



In Einstein-Maxwell (EM) theory the Kerr-Newman (KN) solutions represent stationary asymptotically flat black 
holes, characterized uniquely by their global charges: their mass M, their angular momentum J, their electric charge 
^ I Q, and their magnetic charge P P, @]. Following Wheeler this uniqueness theorem of EM theory is often expressed 
I as "EM black holes have no hair" . 

^ = In many unified theories, including string theory, dilatons appear. When a dilaton is coupled to EM theory, this 
has profound consequences for the black hole solutions. Not only do charged static Einstein-Maxwell-dilaton (EMD) 
black hole solutions exist for arbitrarily small horizon size Q , but also the staticity theorem of EM theory 4] does not 
^ ' generalize to EMD theory for arbitrary dilaton coupling constant 7: at the Kaluza-Klcin value 7kk = \/3 stationary 
^ , non-static black holes appear, whose horizon is non-rotating ^, and beyond 7kk even counterrotating black holes 
, arise, whose horizon angular velocity and global angular momentum have opposite sign 

The EM uniqueness theorem, on the other hand, does not readily generalize to theories with non-Abelian gauge 
fields coupled to gravity The hairy black hole solutions of SU(2) Einstein- Yang-Mills (EYM) and Einstein- Yang- 
Mills-Higgs (EYMH) theory possess non-trivial magnetic fields outside their regular event horizon and are not uniq uely 
characterized by their mass, their angular momentum, their electric and magnetic charge [1, [l^, [HI, [13, [3 Il4{ . 
' Futhermore, black hole solutions arise, which are static and not spherically symmetric, showing that Israel's theorem 
^> . [I, [1| does not generalize to non-Abelian theories, either [l^, [U . 

;V ■ The coupling to non-Abelian fields not only gives rise to new types of black hole solutions, but also allows for 
' globally regular solutions, not present in EM theory either [3, [3, [iHl- These are stationary solutions with a spatially 
localized energy density of the matter fields and a finite mass, and are referred to as solitons when they are stable, 
5-j , and sphalerons when they possess unstable modes. The known globally regular solutions of EYM theory represent 
sphalerons, whereas the globally regular magnetic monopoles of EYMH theory are solitons, whose topological charge is 
proportional to their magnetic charge [l^ . Besides magnetic monopoles EYMH theory contains a plethora of further 
globally regular solutions, representing for instance monopole-antimonopole pairs, chains, and vortex ring solutions 

It is an interesting question whether such globally regular solutions can be endowed with rotation, like their black 
hole counterparts can. When EYM black holes start to rotate, the time component of their gauge potential is excited, 
as expected. Surprisingly, however, not only a magnetic moment is induced by the rotation but also an electric charge 
[T^ [T8!| , and this seems to preclude the existence of globally regular rotating EYM sphalerons [l^, [l^l . 

Globally regular EYMH solutions, on the other hand, can carry electric charge, and the presence of a time component 
of the gauge potential renders the solutions stationary. Together the electric and magnetic fields then give rise to 
an angular momentum density, except in the spherically symmetric case [1^ [2l|, [2^ . Still, globally regular EYMH 
solutions with a non-vanishing global magnetic charge cannot rotate: their angular momentum vanishes [20, 21, 22]. 
But globally regular EYMH solutions with no global magnetic charge do possess a finite angular momentum. In fact, 
it is proportional to their electric charge (20| , giving rise to a quantization condition for the angular momentum, 

J = nQ(l-e), P = ne, e = 1 [1 - (-1)"] , (1) 
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where m and n are two integers, characterizing the EYMH solutions ■ 

Here we derive a mass formula for the stationary globally regular EYMH solutions in the presence of a dilaton. 
Then we address the dependence of the global charges and of the effective action of these solutions on the gravitational 
coupling strength. For a given type of solution, typically two branches of solutions arise, which bifurcate at a maximal 
value of the coupling, a^ax- For static solutions, the mass M exhibits a "spike" at c^max [§,[13) since there the two 
branches must possess the same tangent w.r.t. a [23l |. When stationary and rotating solutions are considered, in 
contrast, the mass branches may exhibit a "loop", when considered as a function of a ^2^. Here we show that for 
stationary and rotating solutions it is the effective action S'^^ which may only exhibit a "spike" in the vicinity of the 
maximal value of the gravitational coupling constant. We illustrate this qualitative different behaviour of the mass 
M and the effective action S'^^ for several sets of numerically constructed stationary and rotating solutions. 

Turning to black holes again, we recall, that EM black holes satisfy the laws of black hole mechanics [l^l and the 
Smarr mass formula [25| 

M = 2TS + 2nj + 4,,iQ + ij,^,^P , (2) 

where T represents the temperature of the black holes and S their entropy, f2 denotes their horizon angular velocity, 
and V'ei and V'mag represent their horizon electric and magnetic potential, respectively. 
In the presence of a dilaton an equivalent mass formula for EMD black holes is 

M = 2TS + 2nj + — + 2iP^iQ , (3) 

7 

with dilaton charge D and dilaton coupling constant 7. Interestingly, this second form of the mass formula also holds 
for the known non-Abelian black hole solutions of Einstein- Yang-Mills-dilaton (EYMD) theory [l^ . 

Here we address stationary black holes of EYMH and Einstein- Yang- Mills-Higgs-dilaton (EYMHD) theory [2^ . For 
these black holes the zeroth law of black hole mechanics holds pj| , as well as a generalized first law . We derive a 
mass formula for EYMHD black holes, based on the asymptotic expansion of the metric and the matter fields. The 
analytical mass formula represents a good criterion for the quality of numerically constructed EYMHD black hole 
solutions, also presented. 

In section II we recall the SU(2) EYMHD action and the equations of motion. We discuss the stationary ansatz 
for the metric, the gauge potential, the Higgs field and the dilaton field, and we present the boundary conditions 
for globally regular and black hole solutions. In section HI wc address the physical properties of the solutions. We 
present the asymptotic expansion at infinity and the expansion at the horizon, needed to obtain the global charges 
and the horizon properties of the solutions. We evaluate the mass, the angular momentum and the effective action of 
the globally regular solutions in section IV, and discuss the dependence of these quantities on the coupling constant 
a. We illustrate these results for a set of numerically constructed solutions. We then derive the mass formula for 
the stationary black hole solutions in section V, presenting also numerical results. In section VI we present our 
conclusions. 



II. EINSTEIN- YANG-MILLS-HIGGS-DILATON SOLUTIONS 



After recalling the SU(2) EYMHD action and the general set of equations of motion, we discuss the ansatz for the 
stationary non-Abelian globally regular and black hole solutions. The ansatz for the metric represents the stationary 
axially symmetric Lewis-Papapetrou metric in isotropic coordinates. The ansatz for the gauge potential and the Higgs 
field includes two integers, m and n, related to the polar and azimuthal angles. For monopole-antimonopole chains 
the integer m counts the total number of poles on the symmetry axis, while the integer n gives the magnitude of the 
magnetic charge of each pole. As implied by the boundary conditions, the stationary axially symmetric solutions are 
asymptotically flat, and the black hole solutions possess a regular event horizon. 

A. SU(2) EYMHD Action 

We consider the SU(2) Einstein- Yang-Mills-dilaton action 
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where R is the scalar curvature, and the matter Lagrangian Lj\/ is given by 



1 o.,„ 1„ A 



d^^d^^ - -e2«*Tr(F^,F'^'^) - -Tr - -e-2«*Tr ($2 _ y^V , (5) 

A Z 4 o 

with dilaton field \I>, gauge field strength tensor = dfj,A^ — d^A^^ + «e [A^, A^], gauge field A^ — A'^Ta/2, Higgs 
field in the adjoint representation $ = r^^", gauge covariant derivative = + ie [A^, ■ ], and Newton's constant 
G, dilaton coupling constant k, Yang-Mills coupling constant e, Higgs self-coupling constant A, and Higgs vacuum 
expectation value v. 

The nonzero vacuum expectation value of the Higgs field breaks the non-Abelian SU(2) gauge symmetry to the 
Abelian U(l) symmetry. The particle spectrum of the theory then consists of a massless photon, two massive vector 
bosons of mass Mw = ev, and a massive Higgs field Mh — \/2J.v. In the limit A = the Higgs field also becomes 
massless. The dilaton is massless as well. 

Including a boundary term [28| . variation of the action with respect to the metric and the matter fields leads, 
respectively, to the Einstein equations 

G^i, ~ — -g^iuR — SttGT^i^ (6) 

with stress-energy tensor 



-f ^Tr (d^'^D,'^ - ]^g^,D^<^D"<^\ - ^g^,e-2«*Tr ($2 - v^)\ (7) 



and the matter field equations, 

1 
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□ * = Ke^-^Tr [F^uF^") ~ ^^e-^-^^Tr ($2 _ y'^f , (9) 



where □^1' = ^.^'"^ , and 



D^,D^'<^ = Ae^^^^Tr ($2 _ -y2^ $ . (10) 



B. Stationary Axially Symmetric Ansatz 

The system of partial differential equations, Eq. ([S]), Eq. ([5]), Eq. ([5]), and Eq. (fTU]) is highly non-linear and 
complicated. In order to generate solutions to these equations, one profits from the use of symmetries, simplifying 
the equations. 

Here we consider solutions, which are both stationary and axially symmetric. We therefore impose on the spacetime 
the presence of two commuting Killing vector fields, ^ (asymptotically timelike) and 77 (asymptotically spacelike). Since 
the Killing vector fields commute, we may adopt a system of adapted coordinates, say {t, r, 6, Lp}, such that 

i = dt, ii^d^ . (11) 

In these coordinates the metric is independent of t and Lp. We also assume that the symmetry axis of the spacetime, 
the set of points where 77 = 0, is regular, and satisfies the elementary flatness condition 

^^ = 1, X^r^, . (12) 

Apart from the symmetry requirement on the metric {C^g = C^g = 0, i.e., (7^^ = g^^{r^9y), we impose that the 
matter fields are also symmetric under the spacetime transformations generated by ^ and 77. 
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This implies for the dilaton field 

r^^- = Z:^* = , (13) 

so 5* depends on r and 9 only. Introducing two compensating su(2)-valued functions and W^, the concept of 
generalised symmetry [lO, [1^ requires for the Higgs field 

= ie[$, W^] , = ie[$, W^] , (14) 

and for the gauge potential A = A^dx'^, 

(£^A)^ = , 

(C^A)^ = D^W^ , (15) 

where and satisfy 

C^Wr, - A, We + [W^, Wr,]=Q. (16) 

Performing a gauge transformation to set — 0, leaves $, A and independent of 

By virtue of the Frobenius condition and the circularity theorem, the metric can then be written in the Lewis- 
Papapetrou form, which in isotropic coordinates reads 

/ / L f J 

where f,m,l and w are functions of r and 9 only. 

The z-axis represents the symmetry axis. The regularity condition along the z-axis Eq. (|12p requires 

m\g=o.7v = ^|9=0,7r • (18) 

The event horizon of stationary black hole solutions resides at a surface of constant radial coordinate, r = rn, and 
is characterized by the condition /(rn, 9) — ^1^- The Killing vector field 

X-e+— ?7, (19) 

is orthogonal to and null on the horizon [23 |. The ergosphere, defined as the region in which ^^^^ is positive, is 
bounded by the event horizon and by the surface where 

- / + sin^ 9ju}^ = . (20) 



For the gauge fields we employ a generalized ansatz [Ij, [13|, [2^ , which trivially fulfils both the symmetry constraints 
Eq. and Eq. pB)) and the circularity conditions, 

C(n,m) {n,m) \ / „ \ (n) 

Bi^^-^+B2^^^\dt + A^{dip-'^dt)+ l-Ldr+il-H2)d9j (21) 

C(n,m) (n,m) \ 

i^3^^^ + (l-i?4)^^^ 1 , (22) 

and the appropriate Ansatz for the Higgs field is then given by [l^, [S^l 

$ = ^; (<i>ir^") + $2r^"'™^) , (23) 

where n and m are integers. The symbols Tr"'™'', t-^"''"-' and r,^"'' denote the dot products of the Cartesian vector of 
Pauli matrices, f = (r^,, Ty, r^), with the spatial unit vectors 

g(n,m) „ (sin(TO0) cos(ni^), sin(m0) sin(n(p), cos(m0)) , 
.(n,m) _ (cos(mf?) cos(niy9), cos(m0) sin(n(^), — sin(TO0)) , 
e(") = (-sin(n(^),cos(n^),0) , (24) 
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respectively. Like the dilaton field function the gauge field functions Bi and Hi and the Higgs field functions 
depend only on the coordinates r and 9. 

The ansatz is form-invariant under Abelian gauge transformations U [lol . [30l | 



U = expi^-T^">Tir,e)j . (25) 

With respect to this residual gauge degree of freedom we choose the gauge fixing condition rdrHi — dgH2 = 0. For 
the gauge field ansatz, Eqs. ([2T|) - ((22|) . the compensating matrix Wr, is given by 

W,^n^. (26) 

C. Dimensionless Quantities 

Let us now introduce the dimensionless quantities, beginning with the dimensionless coupling constants a, /3 and 7 



2^2 



The dimensionless coordinate x is given by 



A = e^P^ , K = 7 . (27) 



r = X , (28) 

ea 



the dimensionless electric gauge field functions Bi and B2 are 



Si^^Si, B2^^B2, (29) 
V47rG' V47rG 



and the dimensionless dilaton function tp is 



Introducing these dimensionless quantities into the EOMs, the resulting equations depend only on the parameters a, 
/3, and 7. Note, that in the limit 7^0 the dilaton decouples and the equations of EYMH theory are obtained. 

D. Boundary Conditions 

Boundary conditions at infinity 

To obtain asymptotically flat solutions, we impose on the metric functions the boundary conditions at infinity 

f\x—c>o — '^\x—oc — ^la; — 00 1 ; ^la:— cxd — • (^-^) 

For the dilaton function we choose 

V'U=oo = , (32) 

since any finite value of the dilaton field at infinity can always be transformed to zero via -0 ip—iplpo) , x — s- xe~^^^°"^ . 
The asymptotic values of the Higgs field functions (f>i are 

$lU=oo = l, $2U=oo=0. (33) 

We further impose, that the two electric gauge field functions Bi satisfy 

Bl\x=oo=y. S2U=oo=0, (34) 

where the asymptotic value v is restricted to < < 1, and that the magnetic gauge field functions Hi satisfy 

-H'lU=oo=0, i/2U=oo = l~TO, (35) 
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cos9 — cos(m6) 1 — cos(m6') 

Hslx^oo = r-7 m odd , iiaU^oc- = r-7 m even , (36) 

sm tl sm w 

-f^4U=cx) = 1 ' ^ r,^ ■ (37) 

sma 

Boundary conditions at the origin 

To obtain globally regular solutions, we must impose appropriate boundary conditions at the origin. Regularity 
requires for the metric functions the boundary conditions 

dxf\x=o = dxm\x=o = dxl\x=o = , uj\x=o = , (38) 

and for the dilaton function 

a,VU=o = , (39) 

the gauge field functions Hi satisfy 

Hi\x=0 = Hslx^o — , H2\x=Q — H4\x=o — 1 , (40) 

while for even m the gauge field functions Bi and the Higgs functions $i satisfy 

[sin(TO6l)$i + cos(to6')$2] L=o = ' (41) 

dx [cos(m^)$i - sin(m0)$2]L=o = > (42) 

[sm{me)Bi + cos{m6)B2 = O] = , (43) 

dx [cos(m0)5i - sin(m0)52] = , (44) 
whereas for odd m they satisfy Bi\x=Q — ^i\x=o — 0. 

Boundary conditions at the horizon 

The event horizon of stationary black hole solutions resides at a surface of constant radial coordinate, x = xh, and 
is characterized by the condition /(xh, 0) = fl2l|. 

Regularity at the horizon then requires the following boundary conditions for the metric functions 

f\x=xn = m\x=xn = Ax=xn = , Lu\x=xn = t^H = COnst. , (45) 

for the dilaton function 

dxMx=x^ = , (46) 
while the Higgs and the magnetic gauge field functions satisfy 

dx'^i\x=x^ ^ dx'^2\x=x^ ^ , (47) 

Hi\x=xii = , dxH2\x=xii = dxH2,\x=xii = dxH4\x=x^ — , (48) 

with the gauge condition doHi = taken into account [12]. The boundary conditions for the electric gauge field 
functions are obtained from the requirement that for non-Abelian solutions the electrostatic potential is constant at 
the horizon [3] 

§ely = -X^^Mlr=rH • (49) 

Defining the dimensionless electrostatic potential ^oi, 



V^el = ^*el , (50) 

a 
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and the dimensionless horizon angular velocity fi, 



n^^, (51) 



yields the boundary conditions 

Bi\x=x-ti = nftcosmO , B2\x=xh = —ni^smmO . (52) 

Boundary conditions along the symmetry axis 

The boundary conditions along the z-axis {9 = and 9 = n) are determined by the symmetries. For the positive 
z-axis they are given by 

defle^o = dem\e=o = del\e=o = deuj\e=o = , (53) 

dey^le^o = , (54) 

i/i|e=o = i?3|e=o = , deH2\e=o ^ deHile^o ^ , (55) 

B2|e=o = 0, 5eBi|e=o = 0, (56) 

*2|e=o = 0, de<^>i\e=o = . (57) 

The analogous conditions hold on the negative z-axis. We note, that the globally regular solutions are symmetric 
w.r.t. the xy-plane. For the black hole solutions, this symmetry is broken via the boundary conditions of the time 
component of the gauge field. 

In addition, regularity on the z-axis requires condition Eq. (jlSp for the metric functions to be satisfied, and regularity 
of the energy density on the z-axis requires 

H2\e=o ^ Hi\g^o ■ (58) 



III. PROPERTIES OF REGULAR AND BLACK HOLE SOLUTIONS 



We derive the properties of the stationary axially symmetric solutions from the expansions of their metric and 
matter field functions at infinity, at the origin and at the horizon. The expansion at infinity yields the global charges 
of the solutions, the expansion at the horizon yields the horizon properties of the black holes. 
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A. Asymptotic Expansion 

The asymptotic expansion depends on the integers m and n. Here we restrict to odd winding number n, since the 
analysis for the even case seems to be 'prohibitively complicated'. We then obtain for /? = 

Cisin6» C3sin(26') ^/l 
Hi = h U 



where 



Cicose C^cos(2^ 
H2 = {l-m) — - + O 

cos( zO )- co^nie) Cesine f 1 
-H3 = ^-7 1 ^ <^ ^ 

H - (1- ^ cos{s0) _ C3 cos^ 9 + C'iCq sin^ 1 

^ \ sin 9 ) X x"^ " 

^ _^ Q ^ 2Q {11 - J D) - lyC fsin^ 9 + 2C7 cos 9 ^ ^/l 
x 2x^ \x^ 

- uCi sin e {CiQ- vCs cos 6») sin 61 ^ / 1 

-D2 = 5 \- O \ ^ 

x \ x-^ 

_ 2^ ^ 2M' + a^(Q' + P')+C4Cos^ + O f — 
X x^ \x^ 

Cs cos(2g) + [-/x^ + a^(Q^ + p2 _ j,2 _ Q'ly^ ^^^2 g 

x^ 

x^ 



D ^{Q^-P^)-2Cscos9 ^ ^ / 1 

X 2x2 \ X'* 



^ C2 C?sin2 6'- 2C9COS6' ^/l 

$. = ^+(^^±%^i^ + 0f^) , (59) 
X x^ \x-^ ' 



£ = i[l-(-in, P = nE. (60) 

For generic /3 solutions, the expansions remain valid with C2 = Cg = 0. At first sight the power law decay of 
the Higgs field then appears surprising, since /3 ^ renders the Higgs massive and should thus lead to an exponential 
decay. However, this power law decay represents a gauge artifact and can be removed by the gauge transformation 

U = exp(irT^/2) , (61) 

with 

r = -(i-r.). + ^ + ^3ff). (62) 

Performing this gauge transformation leads to $ = WTr"'^' + 0{l/r^) (note that the integer m has been transformed 
away). We further obtain trivial gauge field functions Hi and H2 (up to order 0{l/x^)). 
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B. Global Charges 

The expansion coefficients M, J, Q, P and D correspond to the global charges of the solutions. The dimensionless 
mass M and angular momentum J of the solutions are obtained from the asymptotic expansion of the metric 

M = -i^ lim x'^dxf = A , J ^ TTT 1™ ^^'^ = ~^ ' (63) 

These correspond to the expressions obtained from the respective Komar integrals, as shown in sections IV and V 
for the globally regular and black hole solutions, respectively. Note that M is the mass in units of A-Kv/e, whereas jx 
corresponds to the mass in units of Planck mass. Likewise, the dimensionless electric charge Q and the dimensionless 
magnetic charge P are given by 

— 77 

Q = - Yim x{B^-u) , P = - (1 - {-ID , (64) 
while the dimensionless dilaton charge D is given by 

D = lim x^d^^ . (65) 



C. Expansion at the Horizon 

Expanding the metric and matter field functions at the horizon in powers of 

X 



yields to lowest order 



- 1 (66) 
xn 



H2 = H20 + O{S^) , 
= Hso + 0{5') , 
Hi = Hi0 + O{S^) , 

Bi = n— cos {m9) + 0{S'^) , 

Xh 

B2 = -n—sm{me)+0{6^), 

f = 6^2(1 - S) + 0{S') , 

m = 6^m2{l - 36) + 0{6'^) , 

I = S%{1 - 36) + 0{S^) , 

w = u}h{1 + 6) + 0{6^) , 

= + 0{6^) , 

$1 = $io + 0(<5') , 

$2 = <^2o + 0{6^) . (67) 

The expansion coefficients /2, m2, I2, tpo, Hn, H20, H40, $10, $20 are functions of the variable 0. Among these 
coefficients the following relations hold, 

0=^-2%^, (68) 
m2 /2 



Hii = deH2o 



(69) 
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D. Horizon Properties 

With help of the above expansion we obtain the horizon properties of the SU(2) EYMHD black hole solutions. The 
first quantity of interest is the area of the black hole horizon. The dimensionless area is given by 



Ah = 2^ r dO sine^^^xl , (70) 
Jo h 

and the dimensionless entropy S by 

S-^. (71) 
The surface gravity of the black hole solutions is obtained from 

4 = -^(V,.X.)(V^x'') , (72) 

with Killing vector x, Eq. (|19p . Inserting the expansion at the horizon, Eqs. (|67p. yields the dimensionless surface 
gravity 

^ . (73) 

As seen from Eq. (I68p . Kgg is indeed constant on the horizon, as required by the zeroth law of black hole mechanics. 
The dimensionless temperature T of the black hole is proportional to the surface gravity, 

T=^. (74) 



E. Electric and magnetic charge 

A gauge-invariant definition of the electromagnetic field strength tensor is given by the 't Hooft tensor [l^ 

:F^. = Tr |<l^^,,, - ^^D^^dA^ - + ^eatA^'D^^'DA' , (75) 

where 4 is the normalized Higgs field, |$p = (l/2)Tr$^ = Z]a(*")^ = 1- 
The 't Hooft tensor yields the electric current 

\I^T^^ = -4^j-i , (76) 

and the magnetic current 



= 47rj- , (77) 



where * J- represents the dual field strength tensor. 
The electric charge Q is given by 



Q=- = ^ [ *:Fg^d9dip , (78) 

'S2 



e An 



where the integral is evaluated at spatial infinity. 

To define the magnetic charge, we rewrite the 't Hooft tensor as 



:F^. = d^A, - d,A^, ~ ^Tr {^d^dA} , (79) 

with Afj, = Tr I $ I . Now it follows from Eq. (|77)) that the magnetic current and the topological current k"^ are 
related by 

^ "Trldp^d^d,^} ^-k'' . (80) 
loTre 1^ f- J e 
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For globally regular solutions the integration of the magnetic charge density reduces to a surface integral at spacial 
infinity which yields 



V=-£ . 

e 



For black hole solutions we define the magnetic charge by its value on the horizon plus a volume integral, 

V = Vn+ I i-jm^,n^')dV = Vh+ [ jl,^drdedip , (81) 

where E now denotes an asymptotically flat spacelike hypersurface bounded by the horizon H, dV is the natural 
volume element on S, and n'^ is normal to S with n^n'^ = —1. 

In order to define the horizon magnetic charge we consider the normalized Higgs field at the horizon as a map 
between two two-dimensional spheres, which can be characterised by a topological number, 

^ii = ^j Tr|$dl>Ad$| , (82) 

and obtain Vn = Nn/e. For the evaluation of the magnetic charge we note that the volume integral reduces to a 
surface integral. Its contribution from the horizon cancels exactly the horizon magnetic charge, and the contribution 
from the asymptotic region yields V = en/e. Note that for odd m the horizon magnetic charge is either equal to 
the magnetic charge or to its negative value, depending on how often the Higgs field function <I>i changes sign on the 
symmetry axis. For even m both the magnetic charge and the horizon magnetic charge are zero. 



F. Physical Interpretation of v 

The quantity i> is related to the asymptotic behaviour of the gauge potential Aq, and therefore it is not defined 
in a gauge-invariant way. To find a physical interpretation of v we apply a gauge transformation that leads to an 
asymptotically trivial gauge potential (for even m). Such a gauge transformation is given by 

The transformed gauge potential and Higgs field are found to be 



Bi — v — n—{cos{mO) — 1) 

X 



2e 



+ 



B2 + n— sin(TO' 



X 



2e 



dt 



/ TT \ {n,i/t) 

+A^{dip - -dt) + ( —dx + {1-H2- m)d9 ) ^ — , 



(83) 



with 



A^ = — nsinf 



H3 + 



coa{in9) — 1 



sm( 



2e 



1-Ha- 



sin(TO6 



smt 



Ip 

2e 



(84) 



and 



(85) 



respectively, where now 



= cos{nLp — i't)Tx + sm{nip — i't)Ty , 

^(riji^t) _ _ gjjj(^j^(^ _ l^t)Tx + COs{n(p — h't)Ty 



We observe that in this gauge the fields are explicitly time dependent and rotate in internal space about the 
direction. The quantity u is exactly the rotation frequency. 

In the presence of magnetic charge, i.e. odd to, the transformed gauge potential is singular on the negative z axis. 
However, the physical interpretation of ly does not change. 
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IV. STATIONARY GLOBALLY REGULAR EYMHD SOLUTIONS 

A. Global Charges 

Mass, angular momentum and dilaton charge 

We begin by recalling the general expressions [24] for the global mass 

and the global angular momentum 

J^-^ I R^.un^^v'dV . (87) 

Here E denotes an asymptotically flat spacelike hypersurface, n'^ is normal to E with n^n'^' = —1, and dV is the 
natural volume element on S [23 |. 

Now we express the Ricci tensor in terms of the Yang-Mills, Higgs and dilaton fields, using the Einstein equations, 
the definition of the stress energy tensor and the Lagrangian 

^ R,. = d^^d,^ + 2e2''*Tr(F^"i^,„) - ie2-*Tr(Fp,F''-)g^, 



inG 

+iTrp^$i^,$) + ^ 



Next we replace the third and the last term in Eq. ([88|) via the dilaton equation 

= + 2e2«*Tr(F/^^,„) + iTr(D^$D,<I>) - ^^dx{V^d^^)g^, . (89) 

OTTCt z Zk y — g 

Since ^ and rj are Killing vector fields and since rj is tangential to S, we have 

e'^a^^'-O, T^^'d^^^O, nV5M- = 0, (90) 

and consequently. 



^i?^.n^r = 2e2-*Tr(i^^"F.„)n^r + jTr(i?^ci>i?,ci>)n^r - ^^dxiV^d^^)n^i^ , (91) 



^ = 2e2«*Tr(^^^"F,„)n''ry^ + ^Tr(i?^$i?,$)n^ry'' . (92) 

We now define the dilaton charge V via 

/ -^dxiv^d^mn''^^ dV = -inV . (93) 
Making use of the dilaton charge 27, we obtain for the mass M 

M^A f {e2«*Tr(F^"f^,„)n'^r}rf^+ / {MD^'i>D,^)n^C} dV + —V , (94) 
while the angular momentum is given by 

- -2 ^ {e2«*Tr(i^^"F,„)n^77''} dV - ^ {Tr{D^^D,<P)n^'r^''} dV . (95) 
To evaluate the integrals in Eq. ([M|) and Eq. ([95)1 we use local coordinates {t, r, 9, ip). In these coordinates 

= = (1,0,0,0) , = (0,0,0,1) , dV = ^^drd9dip , (96) 
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and we obtain 



M - —V =Im^-^ J e^^^Tr [Fof.F^''] yf^drdOdLp- J Tr [Do$-D°*] yf^drdOdip , (97) 

J = Ij = 2 I e^-^^Tr ^/^drdedip+\ [ Tr [L>^$D°$] ^/^drdOdip , (98) 

defining the integrals Jm and Xj. 

To evaluate the integrals X^, Eq. ([QT]! and Ij, Eq. ([98]) . we make use of the symmetry relations, Eqs. (fTS]) [20l |. 



where = + ie[A^, ■ ]. The integrals then read 



Df.AaF'^'' ^/^drdedip- / Tr [Do*!?"*] ^/^drdOdip 



(99) 



(100) 



Jj = -2 I e2''*Tr 



Z)^ (A^ - W,,) y/^drd0dip+- / Tr [D^^D'^*] y/^drd0dip . 



(101) 



Adding zero to the above integrals, in the form of the gauge field equation of motion for the zero component, we 
obtain 



Tm^^J Tr 



Df, {Aoe2'"*F°^'y^}J drdOdip + j Tr [ieAo[^, - Do$D°$] ^f^drdOdip , (102) 



Ij = -2 [ Tr [l?^ {(A^ - W^) e^'^^F^^V^}] drd^dy^ 

/ Tr [ie (A^ - Wr,) - D^^D"^] y/^drdOdip . 

2 Je 



(103) 



Making use of the explicit form of the ansatz, exploiting in particular Eq. (fH|) . we see, that for both Tm and Tj the 
second integral vanishes identically, leaving only the first integral to be analyzed further. 

Since the trace of a commutator vanishes, we now replace the derivative by the partial derivative in the 
remaining integrals. 



Je 

I.I = -2 / Tr [a^ {{A^ - Wrj) e2''*F°^y^}] drdOdip 



and employ the divergence theorem. The 0-term vanishes, since \/—g vanishes at = and 6 
vanishes, since the integrands a± — Q and p — 2tt coincide, thus we are left with 



Tm = 4/ Tr [Aoe2''*FO'-V^ dOdp , 



(104) 
(105) 

vr, and the (^-term 
(106) 
(107) 



Since the integrands vanish at the origin, the only contributions to Xm and Xj come from infinity. At infinity the 
asymptotic expansion yields to lowest order 



F^^'V^ = Qs'mt 



2e 



oil) 



Aa = D- 



(n.m) 

2e 



A., 



-n sin 6 



-oil), 

cos(et^) — cos(m6') Tr sm(m&) Tg 



sm( 



2e 



sin 9 2e 



oil) 



(108) 
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where 



V4^ 

The integrals Xm and Xj are then given by 



^ = -P=^, e = i(l-(-ir). (109) 



Xm = , Xj = (1 - e) , (110) 



yielding for the mass M and the angular momentum J 



M = —V+—^, J^^^(\-e). (Ill) 



Returning to dimensionless variables, and noting that 



we obtain the mass formula 



-M^^M, J-^,^. V^^, (112) 

^ = a^— + 2a2z/Q A/ = — + 2i^g , (113) 

7 7 

and the quantization condition for the angular momentum Eq. ([T]) 

C = a^nQ{l -e) <=^ J = nQ{l - e) . 

B. Effective Action 

To address the dependence of the globally regular solutions on the coupling constant a, we now consider the effective 
action S'^^ . In particular, we explain the qualitatively different dependence of the mass M for static and for stationary 
solutions. This concerns only such types of regular solutions where two branches of solutions exist. 

For static solutions, the mass M exhibits a "spike" at the maximal value of the coupling Ofmax, where the branches 
merge and end [1, [111, [l3|. The tangent of the mass w.r.t. a must be the same for both branches at amax [23*1. In 
contrast, for stationary solutions, the mass M exhibits a "loop" in the vicinity of the maximal value of the coupling 
Q^max [22l |. Here the tangent of the mass w.r.t. a diverges at amax- The loop is associated with a critical value of a, 
where the two mass branches cross. 



EfTective action and mass 

Let us begin by defining the effective action , 



S""""^ I [ + Lm] ^gSx , (114) 



with the gravitational effective Lagrangian 



IGttG IGttG 
and the matter Lagrangian Lj\/ Eq. ^ 



l6nG 



^ ' ^R-'-^'] (115) 



2 2 4 8 



The divergence term 9^ in the gravitational effective Lagrangian ensures that the varational principle of the effective 
action S'^^ w.r.t. the functions of the ansatz yields the proper set of field equations. For our particular ansatz of the 
metric is given by 

A^" = Vlsm9 (0,r'^dr,dg,0) . (116) 

^ ' ml 
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Reexpressing the curvature scalar R via the Einstein equations, 

R 



SttG ^M*5^* + ^Trp^$i^^$) + ^e-2«*Tr($2_^;2)2 ^ 



then leads to the effective action 



IGTrG 



-9 



z 8 



(117) 



(118) 



Analogously to the derivation of the mass formula, we next replace the two matter terms in via the equation 
of motion of the dilaton field and obtain 



5' 



off 



rd0dip 



(119) 



Since in the local coordinates the second integral agrees with the integral for the dilaton charge Eq. ((93|) and since 
the 0-term in the first integral vanishes, the effective action becomes 



off 



IGttG 



A^'dOdif 



AttV 
2k 



The remaining integral is evaluated with help of the asymptotic expansion of the metric functions, leading to 



= 



which can be rewritten via the mass formula for the regular solutions Eq. (jllip 



^cff 



-M 



Defining finally the dimensionless effective action 5°^ 



S = —a S 



we obtain 



or equivalently 



5" 



off 



(M - i/Q) 



D 



5^=^ = - ( - + j/C 

.7 



(120) 

(121) 

(122) 

(123) 
(124) 

(125) 



Note, that Eq. (|124p for the effective action remains true when the dilaton decouples, i.e. Eq. (|124p also holds for 
EYMH solutions. 



Dependence of the effective action S'^^ and the mass M on the coupling constant a 

To address the dependence of the effective action 5'°'^ on the coupling constant a wc first make the a-dependence 
more explicit. To this end, we express all quantities in the effective action in dimensionless quantities. Denoting the 
by the divergence term corrected dimensionless curvature scalar R (i.e. i? — > i?), the dimensionless matter Lagrangian 
Lm, and the dimensionless determinant of the metric —g, and we obtain the dimensionless effective action. 



'Cff 



1 



Ry/^d^x+ / LM\^d^x 



(126) 



We now take the derivative of w.r.t. a, taking into account that the metric and matter functions, abbreviated 
by Xi, implicitly also depend on a. The derivative has thus two terms 



dS 



cff 



da 



1 

An 



1 f f diRV^ 
167ra2 J I Qx. 



Aa 



ML 



M\ 



-9) \ dX,_^,^ 



da 



(127) 



equations of motion 
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where the second term vanishes for solutions of the equations of motion, and we are left with 

dS''^ 1 



da Sira^ 



Ryf^d-'x . (128) 



From Eq. (|128p we conclude, that the effective action S'^^ must exhibit a "spike" at the maximal value of the 
coupling ttmax, whcrc the two branches of solutions merge and end, since the tangent w.r.t. a must be the same for 
both branches at amax- This conclusion holds for stationary solutions, as well as for static solutions. 

Let us now address the mass M, which is related to the effective action via Eq. (|124p . i.e. 

M = -S"=ff + i^Q . (129) 
Considering the derivative of the mass M w.r.t. a, keeping v fixed, we obtain 

dM dS'^^ dQ 

— = — :^ + '^-r ■ (130) 

da da da 

Clearly, the crucial difference between the tangent of the mass and the tangent of the effective action resides in the 
last term, containing the derivative of the electric charge Q w.r.t. a. It is this term which allows different tangents 
on both branches at amax- In fact, for stationary (non-static) solutions, we observe that this last term diverges at 
ttmax, yielding a divergent tangent also for the mass, as required for a "loop" associated with both mass branches in 
the vicinity of amax- For static solutions, on the other hand, the mass M always exhibits a "spike" at the maximal 
value of the coupling, amax, since the electric charge Q vanishes. 

Note however, that we assumed so far that the value of v is fixed by a boundary condition, Bi\^ — i^. This is in 
contrast to the case where the electric charge is kept fixed by the boundary condition at infinity, 

x^d.,Bi\^ = Q. (131) 

In this case the quantity v is allowed to vary and its value is adjusted by the numerical procedure. More formally, in 
the variation of the effective action a boundary term shows up, which evaluates to QSv. Therefore the field equations 
are obtained from the modified effective action, 

5 (5'=*-Qi/)|Q ==0 (132) 
As a consequence, if a is varied for fixed electric charge Q, 

' S''^ ~Qv)\^^^-^ [ Rs/^d'^x , (133) 



da "3 STra'^ 

since the variation of the modified effective action with respect to the fields vanishes. On the other hand, since 
5<=*^ -Qu = -M, we find 



dM 1 



da %-Ka^ 

Thus for fixed electric charge it is the mass that exhibits a spike. 



Ry/^d-'x . (134) 



C. Numerical results 



We solve the set of thirteen coupled non-linear elliptic partial differential equations numerically J33j, subject to the 
above boundary conditions, requiring the solutions to be regular at the origin. We employ compactified dimensionless 
coordinates, x = x/(l + x). The numerical calculations, based on the Newton-Raphson method, are performed with 
help of the program FIDISOL [33| . The equations are discretized on a non-equidistant grid in x and 6. Typical grids 
used have sizes 100 x 20, covering the integration region < a; < 1 and < 9 < tt/2. (See [H, H^l and for further 
details on the numerical procedure.) 

For given coupling constants a, f3 and 7, the stationary globally regular solutions then depend on the parameter 
v, specifying the time component of the gauge potential at infinity, and on the integers m and n. (In principle, the 
solutions can further depend on the node number of the gauge potential functions fc, labelling the radial excitations. 
However, we here focus on the lowest mass solutions.) 
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We now illustrate the above relations for the mass and the effective action with numerical results presented in 
Figs. 1 and 2. We first consider stationary gravitating dyon solutions with m = 1. Dyons with magnetic charge n = 1 
are spherically symmetric [3ll [ , dyons with higher magnetic charge are axially symmetric [22j . 

In Fig. la we exhibit the mass M of dyons with magnetic charge n — 1, 2 and 3 at a fixed value of ly versus the 
coupling constant a (at /3 = 7 = 0). In each case, a first branch of gravitating dyons emerges from the corresponding 
flat space solution at a = and extends up to a maximal value of the coupling constant, amax, beyond which no dyon 
solutions exist. For the n = 1 dyons we observe a second branch of solutions in the vicinity of the maximal value 
of a. This second branch ends at a critical value of a, where the branch of non-Abelian solutions merges with the 
corresponding branch of extremal Reissner-Nordstrom solutions [sl, |3l|. For n > 1 dyon solutions, numerical accuracy 
does not allow us to discern the existence of two branches. 




— ' — ' — ' — ' — ■ — ' — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — I I ' ' ' ' ' ' ' J 

0.5 1 1.5 0.4 0.8 1.2 1.6 

(c) a (d) a 



Fig. 1 The mass M (a), the scaled mass aM (b), the electric charge Q (c), and the effective action S° (d) are shown versus 
the coupling constant a for dyon solutions with m — 1, n = 1,2,3 and for electrically charged monopole-antimonopole 
resp. vortex ring solutions with m = 2, n — 1,2,3 at ly — 0.32, /3 = 7 = (first branch: solid, second branch; dot-dashed). 



Besides dyons, we also exhibit in Fig. la the mass of electrically charged monopole-antimonopole pair resp. vortex 
ring solutions, which have m = 2 and n = 1, 2 and 3. For these solutions always two branches of solutions exist. 
Again, the first branch emerges from the respective flat space solution, and extends up to a maximal value of a, where 
it merges with the second branch. But the second branch now extends back to a — > 0. The mass diverges on the 
second branch in the limit a 0. But considering the scaled mass aM instead, exhibited in Fig. lb, one realizes, 
that in the limit a — > a globally regular EYM solution ^Td, ^3d\ is reached (after rescaling) [l3| . Clearly, the electric 
charge Q, exhibited in Fig. Ic, also tends to zero on the second branch in the limit a ^ 0, yielding non-rotating 
limiting EYM solutions, in agreement with previous results on globally regular EYM solutions |19l . [20| . 

We exhibit the effective action for the same set of solutions in Fig. Id. As predicted above, the effective action 
exhibits a "spike" close to the maximal value of the coupling constant a, whenever two branches of solutions are 
present. The mass, in contrast, exhibits a "loop" for these stationary non-static solutions close to Ofmax- At ctmax the 
tangent of the mass diverges, since the tangent of the electric charge diverges there. If we consider branches of static 
solutions instead, the mass exhibits a "spike" close to amax H, [l3- 

While in the solutions of Fig. 1 the dilaton is decoupled since 7 = 0, we consider in Fig. 2 the dependence of the 
solutions on the dilaton coupling constant 7. In Fig. 2a we exhibit the mass and the effective action of electrically 
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charged monopole-antimonopole pair solutions {m = 2, n = 1, v = 0.32, /5 = 0) versus the coupling constant a at 
fixed dilaton coupling constant 7 = 0.4. As predicted, we observe a "spike" for the effective action and a "loop" for 
the mass. Note, that the mass of these EYMHD solutions does not diverge on the second branch in the limit a 0, 
since a YMHD solution is approached. We note, that the effect of a dilaton on globally regular solutions is very 
similar to the effect of gravity ^] . 




Fig. 2 The mass AI and the effective action S'^ of electrically charged monopole-antimonopole solutions with m = 2, n = 1, 
u — 0.32, /3 = are shown versus the coupling constant a for dilaton coupling constant 7 = 0.4 (a), and versus the dilaton 

coupling constant 7 for coupling constant a = 0.15 and 0.40 (b). 



The 7-dependence of the solutions (at fixed a) is illustrated in Fig. 2b for electrically charged monopole- 
antimonopole pair solutions with m = 2, n = 1, v = 0.32 (/3 — 0) for a = 0.15 and 0.40. Again we note, that 
the effective action exhibits a "spike" while the mass exhibits a "loop" . That this must be the case, can be shown 
by an argument analogous to the one employed above for the a-dependence. Note, that the mass of these EYMHD 
solutions does not diverge on the second branch in the limit 7 — > 0, since EYMH solutions are approached in this 
limit. 

Wc illustrate the globally regular solutions with an example in Fig. 3. We here exhibit the energy density of the 
matter fields e 

for a monopole-antimonopole pair solution (m = 2, n = 1) carrying electric charge and angular momentum = 0.32, 
a = 0.3, /? = 7 = 0). The maxima of the energy density are associated with the location of the magnetic poles on the 
symmetry axis. 



V. ROTATING EYMHD BLACK HOLES 



A. Non-Abelian Mass Formula 



We now derive the mass formula for stationary axially symmetric EYMHD black hole solutions 

M = 2fS + 2VIJ + — + 2vQ + 2Vm(9(1 - e), (136) 

7 

where T and S are conveniently scaled dimensionless temperature and entropy, respectively. The derivation is analo- 
gous to the derivation of the expressions for the mass and the angular momentum of the globally regular solutions. 
Again we begin by recalling the general expressions [2l| for the global mass 

M = + ^ ^ R^^uU^CdV , (137) 

and the global angular momentum 

^ = Jh - ^ R^,,n^'ll''dV , (138) 
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(a) 

e = 0.2500 



e = 0.3412 



e = 0.3700 



(b) 



Fig. 3 (a) The energy density of the matter fields for an electically charged monopole-antimonopole pair solution (m = 2, 
n = 1, 1/ = 0.32, a = 0.3, /3 = 7 = 0). (b) Also shown are surfaces of constant energy density. 



where S now denotes an asymptotically flat spacelike hypersurface bounded by the horizon H, and the horizon mass 
Mh [24] and the horizon angular momentum J7h are given by 



Ma 



8ttG 



2 m 



IGttG 7h 2 



(139) 
(140) 



Substituting the horizon mass A^h hi Eq. (jl37p and eliminating the horizon angular momentum J7h yields for the 
global mass A4 



M = 2TS + 2—J + 2 



1 



8ttG 



Rpun^'^dV 



1 



(141) 



Following now the same steps as for the globally regular solutions, we obtain analogously to Eq. ([57)1 and Eq. 



M - 2TS - 2— J - —V =1 = 



g2.*Tr 



- / Tr 



y/—gdrd6d(p 
-gdrdOdip , 



(142) 
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defining the integral X. 

We next note, that formally, we can express the integral I in terms of the integrals Eq. ([57)1 and Tj Eq. 
keeping in mind, that S is here bounded by the horizon. Thus 



(143) 



To evaluate X, we proceed again analogously to the globally regular case, making use of the relations obtained for 
Tm and Tj. Analogously to Eq. (|106p and Eq. (|107p we arc then left with 



X = 4 / Tr 



dBdkp , 



(144) 



where the integrand must be evaluated at the horizon and at infinity. Since the electrostatic potential is constant at 
the horizon (see Eqs. (liQ]) ) and 



A,. 



(145) 



the integrand vanishes at the horizon, and the only contribution to X comes from infinity. 

At infinity both expressions have been evaluated in Eq. (|110p . Thus the mass formula becomes 



Al-2r5-2^^J-l^P = ^ + 5^(l-e). 



Returning again to dimensionless variables, recalling Eq. ([50]) and Eq. (|112p . and noting that 



(146) 



(147) 



we obtain the mass formula Eq. (jl36p 
or equivalently 



2TS + 217C + a^ — + 2a2^Q + 2o?^^xQ(\ - e) 
1 
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M = 2X5 + 2VLJ + — + 2vQ + 2V'ciQ(l - e) , 

7 

with scaled dimensionless temperature and entropy, TS = TSj c? . 

This mass formula differs from the EMD and EYMD mass formula Eq. ([3]) in two respects. First, the last term 
is present only for magnetically neutral black holes. Second, the fourth term is an additional term, not present for 
EMD and EYMD black holes. It appears for all electrically charged EYMHD black holes, and has the gauge potential 
parameter v entering together with the electric charge. We note, that the first two terms and the last term do not 
appear in the mass formula Eq. ()113|1 for globally regular solutions. Indeed, when the black hole horizon size is taken 
to zero, the first term vanishes, and the second and the last term cancel, leaving the mass formula Eq. (|113p for 
globally regular solutions. 



B. Numerical Results 



The numerical black hole calculations are performed analogously to the calculations of globally regular solutions 
[33| . except that for black hole solutions we employ the compactified dimensionless coordinate x = 1 — (xu/a;), and 
we impose boundary conditions at the regular horizon. 

For given coupling constants a, /3 and 7, the rotating non-Abclian black hole solutions then depend on the horizon 
radius xh, and on the rotational velocity of the horizon fi in addition to the gauge potential parameter v and the 
integers m and n. 

We exhibit in Fig. 4 an example of a dyonic rotating black hole, which has 771 — 1, n = 1, horizon radius xjj = 0.1, 
horizon angular velocity = 0.5, gauge potential parameter v = 0.04 and the coupling constants a = 0.3, /? = 0.1, 
7 = 0.1. We again exhibit the energy density of the matter fields e, Eq. (|135p . 
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£=0.010 
£ = 0.015 
£ = 0.020 




(a) 
e = 0.010 



e = 0.015 



e = 0.020 





(b) 

Fig. 4 (a) The energy density of the matter fields for a dyonic rotating black hole (m = 1, n = 1, a;H = 0.1, Q — 0.5, f = 0.04, 
a — 0.3, l3 — 0.1, 7 — 0.1). (b) Also shown are surfaces of constant energy density. 



VI. CONCLUSIONS 



We have considered non-perturbative globally regular and black hole solutions of EYMHD theory. These stationary 
axially symmetric solutions are asymptotically flat. The solutions are characterized by two integers, m and n, related 
to the number of monopoles and antimonopoles in the solutions, and to the magnetic charge of the monopoles, 
respectively [r?| . The black hole solutions carry non-Abelian hair outside their regular horizon. 

The globally regular solutions do not rotate, when they carry a global magnetic charge. Only solutions with no 
global magnetic charge can possess angular momentum, which is then quantized in terms of the electric charge [13] ■ 
The globally regular solutions satisfy a simple mass formula 



7 



2vQ 



The presence of electric charge enforces stationarity of the solutions, since it gives rise to an angular momentum 
density (except for the spherical n = 1 monopole). 

The effective action of the globally regular solutions can be expressed in terms of the mass and the electric charge 



5"=*^ = - (M - lyQ) = 



D 

- \ —+vC 
7 
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Based on the effective action we have shown, that the mass of stationary solutions can exhibit a "loop" close to the 
maximal value of the coupling constant amaxi whereas the mass of static solutions can only exhibit a "spike" there 
(for given m, n, /?, 7 and v). 

Rotating EYMHD black hole solutions satisfy the zeroth and the first law of black hole mechanics. Here we have 
derived a non-Abelian mass formula for these black holes, which involves their global charges and their horizon 
properties 

M = 2fS + 2nj + — + 2uQ + 2V;eiQ(l - e) . 

7 

This mass formula differs from the EMD and EYMD mass formula Eq. ([3]), since the last term is present only for 
magnetically neutral black holes, and further an additional term is present for all electrically charged solutions, where 
the gauge potential parameter v is entering together with the electric charge. When the black hole horizon size is 
taken to zero, the mass formula for globally regular solutions is recovered. 

Whether the presence of the dilaton also allows for a new uniqueness conjecture for hairy black holes remains to be 
seen. Clearly, when only the mass, the angular momentum and the electric and magnetic charges are considered, the 
black hole solutions are not uniquely determined by these global charges. 

In the numerical calculations we have only began to investigate the large parameter space for the black hole solutions. 
Here further investigations might reveal new phenomena, not encountered previously for non-Abelian black holes. For 
instance, non-Abelian counterexamples to the staticity theorem might arise as well as counterrotating black holes . 
Also, systems of non-Abelian black holes with regular non-degenerate horizons might exist. 
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